Flow properties of driven-diffusive lattice gases: theory and computer simulation 
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We develop n-cluster mean-field theories (1 < n < 4) for calculating the flux and the gap distribu- 
tion in the non-equilibrium steady-states of the Katz-Lebowitz-Spohn model of the driven diffusive 
lattice gas, with attractive and repulsive inter-particle interactions, in both one and two dimensions 
for arbitrary particle densities, temperature as well as the driving field. We compare our theoretical 
results with the corresponding numerical data we have obtained from the computer simulations 
to demonstrate the level of accuracy of our theoretical predictions. We also compare our results 
with those for some other prototype models, notably particle-hopping models of vehicular traffic, 
to demonstrate the novel qualitative features we have observed in the Katz-Lebowitz-Spohn model, 
emphasizing, in particular, the consequences of repulsive inter-particle interactions. 



I. INTRODUCTION 

The driven-diffusive lattice gas models are of cur- 
rent interest in non-equilibrium statistical mechanics 
jl], H H Depending on the nature of the drive, 

these driven-dissipative systems can attain steady-states 
that are far from equilibrium. The simplest driven- 
diffusive lattice gas model that incorporates inter-particle 
interactions is the Katz-Lebowitz-Spohn model || (from 
now onwards referred to as KLS). Some of the particle- 
hopping models of vehicular traffic |(| are closely related 
to some special limits of the KLS model in one dimen- 
sion. Therefore, in order to compare and contrast the 
spatio-temporal organizations and the flow properties of 
the KLS model with those in the particle-hopping models 
of vehicular traffic, we calculate here those properties of 
the KLS model which are important from the perspective 
of vehicular traffic. 

Over the last decade extensive investigations of vehic- 
ular traffic have been made using the so-called particle- 
hopping models which represent each vehicle by a particle 
H 0, ||] . All these traffic models are defined on discrete 
lattices each site of which, in the spirit of the lattice 
gas models, represents a cell that can accommodate at 
most one particle at a time. In almost all the standard 
particle-hopping models of vehicular traffic the only non- 
vanishing inter-particle interaction is the mutual hard- 
core repulsion which is usually implemented through the 
condition of exclusion principle: no two particles are 
allowed to occupy the same lattice site simultaneously. 
Therefore, a comparison of our results on the KLS model 
with the corresponding results for the particle-hopping 
models of vehicular traffic will show the effects of inter- 
particle interactions other than mere hard-core repulsion. 

The flux (per lane) is defined to be the number of 
particles (per lane) crossing a detector site per unit time. 
In the context of vehicular traffic , the most important 
quantity of interest is the so-called fundamental relation 
which depicts the dependence of the flux on the density 
of the vehicles. The number of empty sites in between 
a pair of particles is usually taken as a measure of the 



corresponding distance- headway DH. 

In this paper we theoretically calculate the DH dis- 
tributions and the flux in the steady-states of the KLS 
model, separately for attractive and repulsive inter- 
particle interactions, within the framework of a cluster 
mean-field theory (MFT) j§ [fi| [□], Q |l|, Q which 
has been very successful also in the theoretical treat- 
ment of the particle-hopping models of vehicular traffic 
H |ll| [l(| |l7). We also indicate the level of the accuracy 
of our cluster-MFT results by comparing these with the 
corresponding numerical data obtained from our com- 
puter simulations of the KLS model. 

The organization of this paper is as follows: in section 
II we define the KLS model and some related particle- 
hopping models which are relevant for our discussion in 
the subsequent sections. We summarize in section III the 
methods of the cluster MFT we use for our theoretical 
calculations as well as those of computer simulation. In 
sections IV, V and VI, we present our theoretical results 
for the one-dimensional KLS model (both with attractive 
and repulsive interactions) in the 1-cluster, 2-cluster and 
4-cluster approximations, respectively, together with the 
corresponding numerical data from our computer simu- 
lations. We present our results for the two-dimensional 
KLS model in section VII. We compare and contrast the 
results for the KLS model with the corresponding re- 
sults for the particle-hopping models of vehicular traffic 
in section VIII before summarizing the main results in 
the concluding section IX. 



II. THE MODELS 
A. The KLS model 

Suppose the variable Cj describes the state of occupa- 
tion of the site i (i — 1,2, ...,N) on a discrete lattice; 
Cj is allowed to take one of the only two values, namely, 
Cj = 1 if the site i is occupied by a particle and = 
if it is empty (or, equivalently, occupied by a "hole"). 
The Hamiltonian for the system, in the absence of any 
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external driving field, is given by 



H 



-4J 



<ij> 



(1) 



where the summation on the right hand side is to be 
carried out over all the nearest-neighbor pairs and J 
takes into account the corresponding inter-particle inter- 
actions. 

The KLS model can be recast in the language com- 
monly used in the theory of magnetism by using classical 
Ising spin variables Si = (2cj — 1) where Si = 1 and 
Si = —1 represent the particles and holes, respectively, 
and the corresponding Hamiltonian, in the absence of the 
external drive, is given by 



J 2_j ^i^j 



(2) 



The attractive and repulsive inter-particle interactions, 
captured by J > and J < 0, respectively, in the Hamil- 
tonian (^) correspond to the ferromagnetic and antiferro- 
magnetic interactions in the form (^|) of the Hamiltonian. 

However, throughout the rest of this paper, we shall 
use the particle-hole picture, where the instantaneous 
state (configuration) of the system at time t is completely 
described by ({c}; t). For example, in case of a system of 
length L in dimension d = 1, ({c};t) = (c.\, C2, C£;t). 
Similarly, for the L x x L y square lattice ({c};i) = 
(cn,ci2, ..,Cij, ...,c Lx L y ;t). The average density c of 
the particles is given by c = hm;v->oo,.£,— «x> N/N s — 
limpf^ 00t L^ 00 (^2f s Ci)/N s where N s , the total number 
of available sites is L for a linear chain and L 2 for a 
square lattice of size L x L. Note that, because of the 
conservation of the particles, the density c is conserved 
by the dynamics. 

The dynamics of the system is governed by the well- 
known Kawasaki dynamics: at any non-zero temperature 
T, a randomly chosen nearest-neighbor particle-hole pair 
is exchanged with the probability min[l, e~^^ H+tE ^\ 
where = (fc^T)" 1 (k B being the Boltzmann con- 
stant) and ATL — Tt({c} new ) ~H({c} id) is the difference 
in the energy of the new and old configurations while 
£ = ( — 1, 0, +1) for jumps, respectively, along, transverse 
to, against the direction of the driving field E. Through- 
out this paper we take fee = 1 and express the tempera- 
ture T in the units of J. 

For the KLS model with attractive inter-particle inter- 
actions (J > 0) on a square lattice, there is not only 
an ordered state at all T < T C (E), but the critical tem- 
perature T C (E) increases with E, saturating at a value 
T C (E -> oo) ~ 1AT C (E = 0) where T C (E = 0) is the 
critical temperature of the corresponding Ising model in 
thermodynamic equilibrium jl], g|. On the other hand, 
T C {E) decreases with E when the inter-particle interac- 
tions are repulsive (i.e., J < 0); the ordering is altogether 
destroyed by sufficiently large E. However, there is no or- 
dered structure at any non-zero temperature in the one- 
dimensional KLS model, irrespective of the sign of the 



interaction J. Because of this intrinsic qualitative differ- 
ence in the nature of the ordering in the steady states in 
d = 1 and d = 2 we present the corresponding results in 
separate sections. 



B. The Totally Asymmetric Simple Exclusion 
Process 

In the Totally Asymmetric Simple Exclusion Process 
(from now onwards referred to as TASEP) pq] , initially, 
N classical particles occupy randomly the sites of a one- 
dimensional lattice of length L (> N). One time step of 
the dynamics consists of updating the position of N par- 
ticles picked up in a random-sequential manner; each ran- 
domly chosen particle moves forward, with probability q, 
if the lattice site immediately in front of it is empty. For 
this model, the simple single-site (i.e., 1-cluster) MFT 
p9[ gives the exact flux 

F = qc(l - c) (3) 

for all densities c of the particles. 

C. Comparison between the models 

In the special case E = the KLS model reduces to the 
corresponding standard Ising model in thermodynamic 
equilibrium. Note that in d = 1, in the opposite limit 
E — > oo the hopping against the driving field becomes 
impossible and, moreover, the field E dominates so over- 
whelmingly over ATt at all non-zero temperatures T that, 
in this limit, the one-dimensional KLS model reduces to 
the TASEP, with q = 1, irrespective of the sign of the 
interaction J, provided J remains finite. However, in 
d — 2, in the limit E — > oo the hopping probabilities in 
the directions transverse to the driving field remain finite 
and the ratio of the hopping probabilities in the direction 
of the field and those in the transverse direction diverges 
pO| . For all the non- vanishing finite E, in the limit J = 0, 
the one-dimensional KLS model reduces to the TASEP 
with the hopping probability q = min{l, e _£B / fcBT }. Fi- 
nally, at infinitely high temperatures each particle moves 
completely randomly, independent of each other, with 
equal probability in all directions. 

III. METHODS OF CALCULATION 

In this section we briefly outline the methods of our 
analytical as well as numerical calculations. 

A. Cluster-mean-field theory 

The dynamical cluster MFT has been used successfully 
in the analytical treatments of several non-equilibrium 
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models including, for example, surface-reaction mod- 
els fllOfl and particle-hopping models of vehicular traffic 
Tfq, [I?]]. However, in all the traffic models there is 
no inter-particle interaction except, of course, the hard- 
core repulsion. Moreover, unlike the traffic models, the 
particles in the KLS model can also move against the 
drive. 

In this paper we extend the approach in appropriate 
manner to calculate the flux F as a function of c in the 
KLS model for arbitrary J^0,0<T<ooandi?^0. 

We define a n-cluster (n < N) to be a collection of n 
sites each of which is the nearest-neighbor of at least an- 
other site belonging to the same cluster. For simplicity 
of notation, let us consider d = 1. We denote the prob- 
ability of finding an n-cluster in the state (ci,C2, ...,c„) 
at time t by the symbol P n (ci, C2, Cn,;t). We treat an 
n-cluster exactly and approximate all the n + m-cluster 
probabilities by a product of n-cluster probabilities in a 
manner so as to couple the n-cluster to the rest of the 
system self-consistently (see, for example, || for a peda- 
gogical introduction and the existing literature). 

It is straightforward to see that the state of the 2- 
clustcr (cj,Cj-|_i) at time t + At depends on the state of 
the 4-cluster (rj_i, Tj, Ti+i, Ti+2) at time t so that the 
exact master equation 



dP 2 {ci,c l+ i) 



dt 



E 



P 4 (Ti_l, Cj+l, a,T i+2 )w{Ti-i, Cj+1, C i: T i+2 ) 



n-i,n+2 

P4,(n-i,Ci, c l+ i,r i+2 )w(r l _i, a, Cj + i,Tj +2 ) 



E 

Tj-l)Ti-2 



P4,(Ti-<i,Ci,Ti-l,Ci+l)w(Ti-2,Ci,Ti-l,Ci+l) 
Pi(Ti-2, Tj— 1) Cjj Cj + i)l0(Ti_2,Tj_i, Ci, Ci+i) 



E 



Pt(Cj, T l+2 , C t+ l,Ti +3 )w(c 2 ,T t+2 ,C i+ l,T i+3 ) 

(4) 



Pi(Ci, C i+1 ,Ti +2 ,T i+3 )w(Ci, C l+1 ,T i+2l T l+3 ) 



governing the time evolution of the 2-cluster proba- 
bilities P2(ci,c,-+i) involves the 4-cluster probabilities 
Pi{T~i-i,Ti,Ti + i,Ti + 2) for all those configurations which 
can lead to the 2-cluster configuration (cj, Cj+i; t) under 
consideration; 

i w(t1,T 2 ,T 3 ,T 4 ) = W(Tl,T2,T3,T4 -> Tl,T 3 ,T 2 ,T 4 ) 

\ e PJ[{ T l T 3+T2T l )-(T 1 T2+T 3 T i ,) + (T2-T 3 )E\^ ^ 

are the corresponding transition probabilities. However, 
the the Master equation governing the time evolution of 
the 4-cluster probabilities P^Tj-i, Tj, 73+1, 1^+2) involve 



6-cluster probabilities, and so on. A few concrete exam- 
ples of such exact master equations for n-cluster proba- 
bilities are given in the appendix A. In the spirit of the 
clustcr-mean-field approach, we truncate this hierarchy 
of exact Master equations by expressing, albeit approxi- 
mately, all the n + m-cluster probabilities in terms of the 
n-cluster probabilities. 

According to the definition of DH, the probability for 
a DH of j is given by 



P(j) = PQJ0000...Q1). 



(6) 



j times 



We evaluate the right hand side of the equation (||) in 
the 2-cluster and 4-cluster approximations. 



B. Computer simulation 

In our computer simulations, we begin with initial con- 
ditions where the particles are arranged randomly on a 
lattice of linear size L (i.e., the system is a linear chain 
of size L in d = 1 and square lattice of size L x L in 
d = 2). The system is then allowed to evolve, follow- 
ing the Kawasaki exchange algorithm with the Metropo- 
lis probabilities min[l, e -/3 ( A7i+£B )] mentioned above, for 
time steps of the order of 10 4 , so as to allow it to reach 
its steady-state and, then, for further 10 5 steps during 
which the time-averaged flux is measured. Finally, for 
the same set of values of the parameters, this process is 
repeated for 100 different initial configurations to com- 
pute the configuration-averaged flux. Since we did not 
observe any significant difference between the results for 
L = 10 4 and L = 10 5 in d = 1 and between those for 
L = 50 and L = 100 in d = 2 all the data presented in 
this paper were generated using L — 10 4 in d = 1 and 
L = 50 in d = 2. 



IV. 1-CLUSTER APPROXIMATION IN 
ONE-DIMENSION 



It has been realized for some time (]J that the smallest 
cluster one must consider in a dynamical cluster MFT de- 
pends on the nature of the dynamics. Since the Kawasaki 
dynamics conserves the number of particles, in principle, 
the smallest cluster must consist of at least one pair of 
sites. However, it is also known that for the E — > 00 
limit of the KLS model, which is exactly identical with 
the TASEP (with q = 1), the single-site MFT gives the 
exact result. Therefore, in this short section we not only 
establish explicitly the limitations of the 1-cluster MFT 
at weak E but also demonstrate how the accuracy of the 
1-cluster MFT increases with increasing E in the KLS 
model in d = 1. 

The net flux is obtained from F = Ff — F r where the 
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forward flux (i.e., flux in the direction of E) is given by 
Ff = ^ ^(cj-i, l,0,Cj +2 )io(ci_i, l,0,Cj +2 ) (7) 

Cj_l,Ci+2 

while the reverse flux (i.e., the flux against £?) is given 
by 

F r = ^ A(ci-i,0,l,Ci + 2)w(cj_i,0,l,Ci +2 ). (8) 

Ci_l,C i+2 

In the 1-cluster approximation the 4-cluster probabili- 
ties P±(ci-\, Ci, Cj+i, are approximated by the prod- 
ucts of corresponding 1-cluster probabilities. Therefore, 
utilizing the facts that P%(1) = c and Pi(0) = 1 — c, in 
the 1-cluster approximation, the equations ^ and (j|) 
give, respectively, 

F f = fi(ci-i)c(l-c)P 1 (c if2 Hc i _i,l,0,c i+ 2)(9) 

C;-l,C i + 2 

and 

F r = J2 ^i(c i -i)(l-c)cP 1 (c i+2 )u;(c i _ 1 ,0,l,c i+2 XlO) 

Ci-l,Cj+2 

for the forward and reverse flux. Hence, the net flux is 



(a) 



F = c(l - c) 



(2c^ - 2c + l)(min[l, e^] - min[l, e~H) 



-PE-i 



c(l — c) I min[l, e 



/3(B-4J)i 



i[l,e^- B+4J >] 



- mm[l,e^- E -^\ 



min[l,e^+ 4J >] 



(11) 



The expression ( |11| ) predicts that the net flux F is inde- 
pendent of the sign of the interaction J, at all c and T, 
irrespective of the strength of E; this is certainly not true 
in general, except at very large E (fig.l). A comparison 
between the theoretical prediction ( |TT| ) and the results 
of our computer simulations (fig.l) exposes not only the 
quantitative inaccuracy of the 1-cluster MFT but also its 
failure to account for the qualitative features of F(c) in 
the case of repulsive inter-particle interactions. 

The fact that the 1-cluster MFT works better for 
stronger E is not surprising as the one-dimensional KLS 
model reduces to the TASEP (with q = 1) in the limit 
E — ► 00 and it is well known that 1-cluster MFT gives 
exact result for the TASEP. Thus, fig-[l] establishes the 
existence of strong correlations which are neglected by 
the 1-cluster MFT. 



V. 2-CLUSTER APPROXIMATION IN 
ONE-DIMENSION 
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FIG. 1: Variation of the net flux F with the density c of 
the particles in the steady-state of the one-dimensional KLS 
model with (a) attractive (i.e., ferromagnetic) interaction 
J = 1.0 and (b) repulsive (i.e., anti-ferromagnetic) interac- 
tion J = -1.0, both at T = 0.5|J|. In both (a) and (b) 
the discrete data points have been obtained from our com- 
puter simulations with E = 1.0(»), 2. 0(+), 3.0(D), 4.0(v) and 
100.0(0), respectively. The lines represent the predictions 
of the 1-cluster approximation for E = 1.0 (solid line), 2.0 
(dotted line), 3.0 (dashed line), 4.0 (long dashed) and 100.0 
(dot-dashed), respectively. 



or, more precisely, 

P4,(n-i,Ti,Ti + i,n + 2) = P 2 (n-i \Tj)P 2 (T il T l+ i)P 2 { n+i \r i+2 ), 

(13) 

Next, we parametrize the 2-cluster probabilities as 
follows: 



In the 2-cluster approximation, we approximate the 4- 
cluster probabilities, appearing in the exact Master equa- 
tion for the 2-cluster probabilities, by a product of 2- 
cluster probabilities, i.e., 

Pi(Ti-l,Ti,T i+1 ,Ti +2 ) OC P 2 (T l -l,T l )P 2 (T il T t+1 )P 2 (T i+ l,T l 

(12) 



P 2 (0,0) = l-c-a (14) 

P 2 (l,l)=c-a (15) 

where 

P 2 (l,0) = P 2 (0,l)=a; (16) 
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this parameterization enables us to satisfy the three con- 
straints Hi namely, £c l+1 =o ft(l, Q+i) = ft(l) = c, 
Ec 1+1 =oWc i+ i) = Pi(0) = 1 - c, and ft (1,0) - 
ft(0, 1). Equivalently, these constraints also imply 
that dP 2 (0,0)/dt = dP 2 {l,l)/dt = -dP 2 {0,l)/dt = 
—dP2(l,0)/dt. So, only one of the four equations rep- 
resented by equation (^) can be taken as an indepen- 
dent equation. Thus, the calculation of the four 2-cluster 
probabilities boils down to the calculation of the single 
parameter a. 

Using Mathematics for an automated generation of the 
equations and simplification, we find that, in the steady 
state, a satisfies the quadratic equation 

Ao 2 + Ba + C = 0, (17) 

where 

A = min[l,e^(- iS - 4J )]+min[l,e' 3 ( B - 4J )] 

- min[l, e «- B +* J )] _ mm[l,e^ E+i \ (18) 

B = — ^min[l, e ^- E - iJ ^} + min[l, e^-*^ , (19) 

C = c(l-c){min[l, e^- E - 4J >]+min[l, e^" 4 ^]}, (20) 

and, hence, taking the physical solution that allows the 
2-cluster probabilities to be between and 1, we get 

a = (-B - ^B 2 - 4AC)/(2A). (21) 

The corresponding form of P2(l,0) for an alternative 
choice of w was derived by Szabo et al.Jl2]|. Both the 
forms ( ^l|) and the equation (13) in reference JlJ] are 
special cases of the general form [pTJ 

a = (1- 2c(! - c)k)/k (22) 

where 

n = 2\l- +E ) I 

[ w(-AH,-E)+w(-AH,+E)y v ' 

w(AH, ±E) being the hopping probabilities against and 
along the field E, respectively. 

Finally, in the 2-cluster MFT, the forward flux and the 
reverse flux are given by 

Ff= J2 P2(Ci-l\r)P 2 (l,0)P2{0\c i+ 2) 

C»_l,Cs-)-2 

xw(ci-i,l,0,Ci + i) (24) 

and 

F r = ft(ci-i|0)ft(0,l)fta|ci+a) 

Ci-l,C i + 2 

xw( Cl _i,0, l,c J+2 ), (25) 

respectively; the net flux is obtained from F = Ff — F r . 

In the next few subsections we compare the predictions 
of this 2-cluster MFT with the corresponding numerical 
data obtained from our extensive computer simulations 
of the KLS model. 



A. Fundamental diagrams for arbitrary E at T > 



It is straightforward to verify that at an infinitely high 
temperature the expressions ( f24|) and (^5|) for Ff and F r 
become identical and, therefore, the net flux vanishes; 
hign temperature not only washes away the effects of J, 
as it is known to do even in equilibrium, but also the 
effects of E making particle movement in both directions 
equally probable. 

In fig. 2 we plot the flux as a function of the particle 
density c, at a fixed non-zero finite temperature T, for 
five different values of E; the agreement between the the- 
oretical prediction and computer simulation is very good 
for both attractive (fig. 2a) as well as repulsive (fig. 2b) 
interactions. 

The F(c) curves are the analogues of the fundamental 
relations for the traffic models. The shape of the curve 
F(c) in the KLS model with attractive (J > 0) is quali- 
tatively similar to those observed in the particle-hopping 
models of vehicular traffic. In sharp contrast, we find 
a qualitatively different shape of the curve F(c) in the 
KLS model with repulsive interactions (J < 0); there is 
a minimum, rather than maximum, at c — 1/2 provided 
the strength of J is comparable to that of E. Moreover, 
for the same T and E, the flux is higher in the case of 
repulsive inter-particle interactions than that in the case 
of attractive inter-particle interactions. Nevertheless, be- 
cause of the particle-hole symmetry, the curves in both 
the figures ||(a) and (b) are symmetric about c = 1/2 
irrespective of the sign of the inter-particle interactions. 

Also note that in both the figures ||(a) and (b) F(c = 
1/2, E) — > 1/4 as E — > 00. This is a consequence of the 
fact that, as stated before, the KLS model reduces to the 
ASEP, with q = 1, in the limit E — > 00 (in both the cases 
of attractive and repulsive inter-particle interactions) as 
long as J and T remain finite. In fact, in all the figures 
0(a), (b) and |(a), (b), the full curve F(c,E -> 00) is 
given by exact expression F(c, E — > 00) = c(l — c). 

Thus, unlike the 1-clustcr MFT, the 2-cluster MFT re- 
produces the qualitative features of the F(c) curves for all 
E and for both attractive as well as repulsive interactions. 
Moreover, the flux predicted by the 2-cluster MFT is also 
in good quantitative agreement with the corresponding 
computer simulation data, except for a narrow range of 
c about c = 1/2. This indicates that the predictions of 
2-cluster MFT, although quite accurate for arbitrary E, 
is not exact (except, of course, E — > 00) for the KLS 
model in d = 1. We shall improve our theory further by 
developing a 4-cluster MFT in the next section. 

In order to get a deeper insight into the dependence of 
the flux F(c, E, T) on c, E, T as well as on the sign of J 
we analyze the results of the 2-cluster MFT in detail in 
a few special limits in the following subsections. 
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FIG. 2: Variation of the net flux F with the density c of the 
particles in the steady-state of the one-dimensional KLS model 
with (a) attractive interaction J = 1.0 and (b) repulsive in- 
teraction J = -1.0, both at T = 0.5| J\. In both (a) and (b) 
the discrete data points, have been obtained from our com- 
puter simulations with E = 1.0(»), 2.0(+), 3.0(D), 4.0(v) and 
100.0(0), respectively. The lines represent the corresponding 
predictions of the 2-cluster approximation. 



B. Flux at T = 

Let us investigate the dependence of the flux 
F(c, E, T = 0) on the driving field E at T = for arbi- 
trary values of c. For repulsive inter-particle interactions, 
the 2-cluster expressions (|l8|)-(|2"o|) for A, B and C reduce 
to the simple forms A = 2, B = —2, C — 2c(l — c) for 
all E < 4|J|, A = 1, B = —2, C = 2c(l - c) at E = 4|J| 
and to the forms A = 0, B = —1, C = c(l — c) for 
all > 4| J|. Substituting the zero-temperature values 
of A, B and C into the equation (|2l|), we find that, at 
T = 0, for all E < A\J\ a = P 2 (1,0) = c for c < 1/2, and 
a = P2(l,0) = 1 — c for c > 1/2; physically, this means 
that when less than half of the sites are occupied by par- 
ticles (holes) , both the nearest-neighbors of each particle 
(hole) are certainly holes (particles) because of the repul- 



sive nature of the inter-particle interaction. On the other 
hand, at T = 0, for all E > 4| J\, a = P 2 (l, 0) = c(l - c) 
for any arbitrary c; this implies that at T = the ef- 
fects of all E > A\J\ is equivalent to those of infinitely 
large £ at all finite non-zero T so that 1-cluster MFT 
becomes exact. Moreover, at T = 0, for _E = 4|J|, 
a = P 2 (l, 0) = 1 - v /l-2c(l-c) for all c. 

Substituting the appropriate expression of a, derived 
above for T = 0, into the limiting form of the net flux 
P(c, E,T — 0), obtained from the 2-cluster expressions 
(e3) and (|2^), we get 



P(c,£,T = 0) = 



F£(c) for E < 4|J|, 
PI(c) for E = 4\J\, 
F^(c) for £ > 4|J|, 



(26) 



for the KLS model with repulsive inter-particle interac- 
tions at T — 0, where 



FZ(c) =6(0.5 -c) 



+0(c-O.5) 



1 -c 
1 - c 



(l-2c) 



(2c -1) 



(27) 



and 



2{2c(l - c) - 1} + 2 v /l-2c(l-c){l - c(l - c)} 



c(l 



P;( C ) = c(l-c); 



(28) 



(29) 



6*(z) being the step function, namely, 6(x) — for all 
x < and = 1 for all x > 0. 

In the special case of half-filling, the particles re- 
main "pinned" to their respective positions by the inter- 
particle interactions J and no E < 4| J\ is strong enough 
to cause any "de-pinning" . Such "switching" of the flux 
from zero to a non-zero value discontinuously, in response 
to the driving field E, has been reported earlier for the 
KLS model on a linear chain jL2| as well as on a square 
lattice 

In sharp contrast to the c-dependence of P<(c) in 
the case of repulsive inter-particle interactions, we have 
P<(c) = for all c when the inter-particle interaction is 
attractive] therefore, in the latter case, 



F(c,E,T = 0) 



F£(c) for E < 4|J|, 
f£{c) for E = A\J\, 
F°{c) for E > 4|J|, 



(30) 



for the KLS model with attractive inter-particle interac- 
tions at T = 0, where 



F<(c)=0, 



{2c(l-c) + l}- v /l + 4c(l-c) 
2c(l - c) 



(31) 



(32) 
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and 



(a) 



F-(c)=c(l-c). 



(33) 



At T = 0, all E > 4| J| is equivalent to E — » oo, irrespec- 
tive of the sign of J and, hence, the corresponding flux 
is c(l — c). 



C. Temperature-dependence of flux 

In this subsection we consider the dependence of F 
on T at a few special values of c. As demonstrated in 
the figures the predictions of the 2-cluster theory 

agree very well with the corresponding numerical data 
obtained from our computer simulations. 

The non-monotonic variation of the flux with temper- 
ature at small and intermediate values of E is an inter- 
esting phenomenon. So long as T is non-zero but much 
smaller than E, it works against J and helps in "de- 
pinning" the particles which can move forward under the 
influence of the driving field E. However, when T be- 
comes much larger than E, then it washes out the effect 
of E allowing particles to move against E as often as 
along E. 



D. DH distribution 

In the 1-cluster approximation, the DH distribution is 
given by 

P lc (j)=c(l- C y. (34) 
However, in the 2-cluster approximation, we write 

ftcO') =Pi (III), fori = (35) 

and 

ftctf) = J P 2 (l|0){P 2 (0|0)y- 1 P 2 (0|l), for j > 1. (36) 
Hence, in the 2-cluster approximation, we get 

P 2c (j) = l--, forj = (37) 
c 

whereas 



PM = 



c(l - c) 



1 - 



1 



j'-i 



, for j > 1. 



(38) 



The variations of the DH distribution with E (for fixed 
c) and with c (for fixed E) , as predicted by (^7|) and j38| ) 
of the 2-cluster approximation, are compared with the 
corresponding computer simulation data in figj^. and 
figj|, respectively. Note that for c = 0.5, in the absence 
of E, the most probable DH is j = or j = 1 depending 
on whether the interaction is attractive (i.e., ferromag- 
netic, in the language of magnetism) or repulsive (i.e., 
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FIG. 3: Variation of the net flux F with the temperature T in 
the steady-state of the one-dimensional KLS model with (a) 
attractive interaction J = 1.0 and (b) repulsive interaction 
J = —1.0, both at c = 0.5. In both (a) and (b) the discrete 
data points, have been obtained from our computer simula- 
tions with E = 1.0(«),2.0(+), 3.0(D), 4.0(v) and 100.0(0), 
respectively. The lines represent the corresponding predic- 
tions of the 2-cluster approximation. 



antifcrromagnetic). This type of spatial organization of 
the particles persists even in the presence of E as long as 
E is much weaker than the strength of the interaction J 
(see fig.^j). However, deviation from this spatial organi- 
zation increases gradually with increasing strength of E. 
In the limit E — > oo, for all finite | J|, the DH distribution 
approaches the exact DH distribution of the TASEP and, 
as expected, is independent of the sign of the interaction 
J. For a given E, which is comparable with the strength 
| J| of the interaction, increasing c leads to more conges- 
tion and, therefore, the probability of having a DH = 
becomes larger for higher c (se fig-H). 
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FIG. 4: Same as in fig.(|3j), except that the density of the 
particles is c = 0.25. 



VI. 4-CLUSTER APPROXIMATION IN 
ONE-DIMENSION 

The exact master equation for the 4-cluster probabil- 
ities (given in the appendix B), as expected on general 
grounds, involve 6-cluster probabilities. In the 4-cluster 
approximation we break up the 6-cluster probabilities in 
terms of the products of the 4-cluster probabilities using 
the prescription 

P6(n-2,Ti-l,Ti,Ti + i,Ti + 2,Ti+3) = Pj(.Tj-2, Tj-1 \ n, Tj+l ) 

xP^n-i, n, n + i, T i+2 ) 

X P 4 ( Tj , Tj+l | T i+2 , T l+3 ) . (39) 

Due to the constraints of conditional probabilities and 
conservation of total probabilities, at 4-cluster level, we 
have 7 independent equations and the same number of 
variables. The solution is obtained by solving a set of 
non-linear equations. We find that although the 4-cluster 
results are, clearly, an improvement over the 2-cluster re- 
sults, the difference in the actual numerical values of the 
flux in the two approximations, for the same set of param- 
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FIG. 5: The distance-headway distribution in the steady-state 
of the one- dimensional KLS model with (a) attractive interac- 
tion J = 1.0 and (b) repulsive interaction J = —1.0, both at 
for c = 1/2 and T = 0.5]J|. In both (a) and (b) the discrete 
data points have been obtained from our computer simula- 
tions with E = 1.0(«),2.0(+), 3.0(D), 4.0(v) and 100.0(0), 
respectively. The lines represent the corresponding predic- 
tions of the 2-cluster approximation. 



eters, is extremely small. The flux obtained in these two 
approximations for a typical set of values of the parame- 
ters are compared in fig|7] (the corresponding differences 
in the DH distributions are too small to be shown in a 
figure) . 

It is worth pointing out that the 4-cluster approxima- 
tion is not only a quantitative improvement over the 2- 
cluster approximation. The 4-cluster approximation can 
account for the forward-backward symmetry breaking, 
an interesting phenomenon |lj, which the 2-cluster 
approximation fails to capture. 
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FIG. 6: The distance-headway distribution in the steady-state 
of the one- dimensional KLS model with (a) attractive interac- 
tion J = 1.0 and (b) repulsive interaction J = —1.0, both at 
for E = 4.0 and T = 0.5|J|. In both (a) and (b) the discrete 
data points have been obtained from our computer simula- 
tions with c = 0.1(«), 0.25(D) and 0.5(v), respectively. The 
lines represent the corresponding predictions of the 2-cluster 
approximation. 
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FIG. 7: Comparison of the predictions of 2-cluster and 4- 
cluster MFTs on the variation of the net flux F with (a) the 
density and (b) temperature in the steady-state of the one- 
dimensional KLS model with repulsive interaction J = —1.0. 
The parameters in (a) are T = 0.5, E — 3.0 while those in 
(b) are c = 0.5, E = 4.0. The dashed and solid lines are the 
theoretical results obtained in the 2-cluster and 4-cluster ap- 
proximations, respectively, whereas the discrete data points 
are the numerical data obtained from our computer simula- 
tions. 



VII. CLUSTER APPROXIMATIONS IN TWO 
DIMENSION 

It is well known jlj] that, usually, the cluster MFT fails 
to account for the properties of the driven-diffusive lattice 
gases below the ordering temperature T C (E). Therefore, 
we confine our discussions in this section to temperatures 
T > T C [E). 

As in d = 1, for the same set of parameters, the 
flux in the two-dimensional KLS model with attractive 
inter-particle interactions (fig.||(a)) is lower than that in 
the same model with repulsive inter-particle interactions 
(figfb)). 

A comparison of the predictions of the cluster MFT 
with the computer simulation data (figj|) establishes 



that, in the case of attractive interactions both the 1- 
cluster and 2-cluster MFT over-estimate the flux, al- 
though the prediction of the 2-cluster theory is closer to 
the computer simulation data. On the other hand, in the 
case of repulsive interactions, the 1-cluster MFT gives an 
under-estimate whereas the 2-cluster MFT provides an 
over-estimate of the flux. The level of accuracy of the 
2-cluster MFT can be estimated from the plots in fig.(||); 
a close inspection, thus, reveals that the 2-cluster MFT 
does not reproduce the dip in the flux around c = 1/2 in 
the case of repulsive inter-particle interactions. 
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FIG. 8: Variation of the net flux F with the density c of 
the particles in the steady-state of the two-dimensional KLS 
model with (a) attractive interaction J — 1.0 and (b) repul- 
sive interaction J = -1.0, both at T = 1.5T C (E = 0). In 
both (a) and (b) the discrete data points have been obtained 
from our computer simulations with E = 1.00) and 4.0(v), 
respectively. In both (a) and (b) the dotted and dashed lines 
represent the predictions of the 1-cluster MFT for E — 1.0 
and 4.0, respectively, while the long-dashed and solid lines 
represent the corresponding predictions of the 2-cluster MFT. 



VIII. COMPARISON WITH THE RESULTS FOR 
OTHER MODELS 

It is well established Q that the 1-cluster MFT result 
(Eq.||) is the exact expression for the flux in TASEP. If the 
random-sequential updating of the TASEP is replaced 
by the parallel updating it becomes identical with the 
Nagel-Schreckenberg (NS) model 0] of vehicular traffic 
with Vmax = 1; where V max is the largest (integer) speed 
allowed for each of the vehicles. 

In the case of the NS model the 1-cluster MFT makes 
an under-estimate of the flux but the 2-cluster MFT 
treatment is adequate to take into account the correla- 
tions introduced purely by the parallel dynamics and, 
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FIG. 9: Comparison of the predictions of 1-cluster and 2- 
cluster MFTs on the net flux F in the steady-state of the 
two-dimensional KLS model with (a) attractive interaction 
J = 1.0 and (b) repulsive interaction J = 1.0, all for E = 
1.0. In both (a) and (b) the discrete data points have been 
obtained from our computer simulations at T = 1.25T C (E = 
0)(o) and T = 1.5T C (E = 0)0), respectively. In both (a) and 
(b) the dotted and solid lines represent the predictions of the 
2-cluster MFT at T = 1.25T C (E = 0) and T = 1.5T C {E = 0), 
respectively. 

therefore, it gives the exact result Jig ]. 

The KLS model can be regarded as an extension of the 
TASEP by incorporating non- vanishing inter-particles in- 
teractions through non-zero J. Our results reported in 
this paper show that neither the 1-cluster MFT nor the 2- 
cluster MFT yield exact flux in the KLS model; although 
the 2-cluster results are accurate to order 10~ 3 it is the 4- 
cluster results that are practically indistinguishable from 
the corresponding computer simulation data. Although 
it is possible that the results of the 4-cluster MFT might 
be exact for the KLS model we refrain from making such 
a claim as we do not have any rigorous proof. 

Note that in the special case 

E = 
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Distance-headway 



sive (anti-ferromagnetic) interactions so long as E is not 
much stronger than \ J\; as E increases the depth of the 
well at c = 1/2 in fig. g decreases and, eventually, for suffi- 
ciently large E the flux exhibits its maximum at c = 1/2. 
Moreover, the location of the maximum in the DH dis- 
tribution depends crucially on the sign of the interaction 
J provided E is not much larger than \J\. 

The time-interval between the arrivals (or departures) 
of the successive particles at a detector site is defined 
to be the corresponding time-headway TH. The exact 
TH distribution for the NS model with V max = 1 has 
been calculated ]17| . However, because of the possibility 
of hopping of the particles against E in the KLS model, 
the analytical calculation of TH distribution is extremely 
difficult and will not be reported here. 

IX. CONCLUSION 



FIG. 10: The discrete data points represent the distance- 
headway distribution in the steady-state of the one- 
dimensional KLS model with repulsive interaction J = —1.0, 
for c = 1/2 at three different values of T = 4 J/ (hip) cor- 
responding, respectively, to p = 0.1 (•), p = 0.25 (□) and 
p — 0.5 (v)> a H f° r E = 0.0. The lines are merely guides to 
the eye connecting the points corresponding to the DH dis- 
tributions in the NS model with Vmax = 1 and p = 0.1 (solid 
line), p = 0.25 (dotted line) and p = 0.5 (dashed line), respec- 
tively. All the results in this figure have been obtained in the 
2-cluster approximation. 



/3J=(l/4)lnp (with. < p < 1), (40) 

we have A = 2q, B = — 2 and C — 2c(l — c) where 
q = 1 — p. In this case, the quadratic equation (|l8| ) 
reduces to the form qa 2 — y + c(l — c) = which was 
derived []15[ directly from the 2-cluster MF treatment of 
the NS model with V max = 1- Interestingly, the 2-cluster 
result is not exact for the KLS model, but it gives exact 
result for the NS model. This is a consequence of the fact 
that the 2-cluster MFT gives, in general exact results for 
the one-dimensional Ising model in equilibrium. 

Since p < 1, the relation ( fi0"| ) maps the NS model, 
with V max = 1, onto an antiferro- magnetic Ising model, 
(or, equivalently, to the KLS model, with repulsive inter- 
particle interactions, in the absence of external drive) so 
that the steady-state of the former is identical to the equi- 
librium state of the latter. Consequently, for all densities 
c, we observe perfect agreement of the DH distributions 
in the NS model and that in the KLS model with E = 0, 
/3J=(l/4)lnp(figjg). 

Since the relation ( pp[ ) maps the NS model (with 
V max — 1) onto the KLS model only for E = 0, this 
mapping cannot relate the properties of the KLS model 
for any non-zero E with those of the NS model. In- 
terestingly, in the NS model with V max = 1 the flux is 
maximum at c = 1/2 for all q. In sharp contrast, the flux 
is minimum at c = 1/2 in the KLS model with repul- 



In this paper we have reported the results of cluster- 
mean-field theoretic treatments of a driven-diffusive lat- 
tice gas model to calculate the flux of the particles under 
the influence of the driving field. Although we have con- 
sidered the standard model, namely, the Katz-Lebowitz- 
Spohn model in this paper, our technique is sufficiently 
general that it can be used, in principle, to calculate the 
flow properties of any other driven-diffusive lattice gas. 

We have shown interesting novel features of the flux- 
density relation of the KLS model with repulsive inter- 
particle interactions. Our investigation has helped in elu- 
cidating the roles of inter-particle interactions J, temper- 
ature T and the driving field E in determining the trend 
of variation of the flux with the density of the particles 
in the driven-diffusive lattice gas models. 
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X. APPENDIX A 

Concrete examples of exact master equations for 
rt-cluster probabilities in the onc-dimcnsional KLS 
model: 



= -min[l,e- /3(B+4J) ]P4(0, 0,1,1) 

+ min[l,e' 3 (- E+4J )]P4(0, 1,0,1) 

+ rmn[l,e^- E+4 ^]P 4 (l ) 0,1,0) 

- min[l,e /3(iS - 4J) ]P4(l,f,0,0), (41) 
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dP 3 (0,0,0) 
dt 



-min[l,e- /JB ]fls(0,0,0,l,0) 
min[l,e-^ +4J )]P 5 (0,0,0,l,l) 



+ min[l,e-Hfls(0,0,l,0,0) 
+ nrfn[l,e^]i%(0,0,l,0,0) 

+ min[l,e' 3 ( B+4J 5]fl5(0,0,l,0,l) 

- min[l, e^Ps (0,1,0,0,0) 

+ min[l, e' 3 (- £+4J )]P5(l, 0, 1, 0, 0) 

- min[l,e' 3(iS - 4J) ]P5(l,l,0,0,0). (42) 

XI. APPENDIX B 

Exact master equations for the 4-cluster probabilities 
in the one-dimensional KLS model: 



dP 4 (ci_i, Ci, c i+ i,c i+2 ) 



ill 



E 



x w(a l ^i,c l+1 ,Ci,c i+2 ) 

— P&{ T i-2, Ci—li Ci, Cj+1, Cj+2, Tj+3) 

x u;(ci_i,Ci,Cj + i,Cj+2) 



E 



Ti-2,T; + 3 



^6(^1-2, Ci-1, C i; Ci+ 2 , Cj+l, Ti+3) 



x u;(ci_i,Ci,Cj+2,Cj+i) 

— P6( r i-2, Ci-1, Ci, Ci + i, Ci+2, Ti+3) 

x u;(ci_i,Ci,c< + i,Cj + 2) 
+ E Pe(Ti-2, Ci, Ci_i, Ci + i, c i+2 , n+3) 

X w(Ci,Ci-i,Ci + i,Ci +2 ) 

— Pe\Ti-2, Ci-i, Ci, Ci+i, Ci + 2, Ti + 3) 

x u;(ci_i,Ci,Ci+i,Cj + 2) 



E 



Pe(cj— 1, Ci, Ci+i, Ti + 3, c.i + 2, Tj+4) 

X w(Cj + i, Ti +3 , Cj+2, Ti +4 ) 

— Pe(Ci-l, Ci, Ci_)_l, Ci+2, Ti +3 , T l+i ) 

X w(c 2+ l,Ci + 2,T l+3 ,T i+i ) 



E 



P6(Ti-3, Ci-1, Ti-2, Ci, Cj+l, Ci+2) 

w(ri_3,Ci_i,T t _ 2 ,C l ) 
P&{n-3-,Ti-2,Ci-l,Ci, Cj+l, Cj+2) 



^(Ti-3, T,-2, Cj-l,Cj) 



(43) 
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